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Cn I Abstract. Precise data for quasi-free photoproduction of 77 mesons off the deuteron have been measured at 

r~^ . the Bonn ELSA accelerator with the combined Crystal Barrel/TAPS detector for incident photon energies 

("^ ' up to 2.5 GeV. The r;-mesons have been detected in coincidence with recoil protons and neutrons. Possible 

^— ( I nuclear effects like Fermi motion and re-scattering can be studied via a comparison of the quasi-free reaction 

T— ( ■ off the bound proton to ry-production off the free proton. No significant effects beyond the folding of the 

IL" ' free cross section with the momentum distribution of the bound protons have been found. These Fermi 

motion effects can be removed by an analysis using the invariant mass of the r;-nucleon pairs reconstructed 

from the final state four-momenta of the particles. The total cross section for quasi-free r;-photoproduction 

' 5^ ' off the neutron reveals even without correction for Fermi motion a pronounced bump-like structure around 

C^ ' 1 GeV of incident photon energy, which is not observed for the proton. This structure is even narrower in 

the invariant mass spectrum of the ?7-neutron pairs. Position and width of the peak in the invariant mass 

spectrum are W « 1665 MeV and FWHM F ~ 25 MeV. The data are compared to the results of different 

models. 

PACS. 13.60.Le Meson production - 14.20.Gk Baryon resonances with S=0 - 25.20.Lj Photoproduction 
reactions 

1 Introduction strong interaction at the few-GeV scale, where perturba- 

tive methods cannot be applied. It plays a similar role for 
The excitation spectrum of the nucleon is one of the most ^^e strong interaction as atomic spectra do for the electro- 
important testing grounds for our understanding of the magnetic mteraction. However, so far, on the theory side 

the only direct connection between baryon properties and 
Correspondence to: B. Krusche, Klingelbergstrasse 82, CH- Quantum Chromodynamics (QCD) has been established 
4056 Basel, Switzerland, e-mail: Bernd.Krusche@unibas.ch "^'^^ ^^e numerical methods of lattice gauge theory. The 
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progress in this field was tremendous during the last few 
years for the ground state properties of hadrons [1] . But 
only very recently also first results for excited states [5] go- 
ing beyond earlier quenched approximations [31I3] became 
available. 

In a more indirect way, experimental observations and 
QCD are connected via QCD inspired quark models. How- 
ever, in spite of their phenomenological successes, the ba- 
sis of these models is still not well anchored. There is nei- 
ther consent about the effective degrees of freedom nor 
about the residual quark - quark interactions (see e.g. 
Ref. [5l[6] for detailed reviews). Apart from the standard, 
non-relativistic quark model with three equivalent valence 
quarks, also models with a quark-diquark structure (see 
e.g. [7]), algebraic approaches 8,, and flux-tube models 
[5] have been proposed, all with different internal degrees 
of freedom, giving rise to different excitation schemes. In 
a radically different picture of hadrons, it was even at- 
tempted to generate all excited states by chirally coupled 
channel dynamics [lU] , leaving only the ground-state mul- 
tiplets as genuine qqq states. 

So far, comparison of the experimentally known exci- 
tation spectrum of the nucleon to model predictions does 
not clearly favor any of the different models, but reveals 
severe difficulties for all of them. The ordering of some of 
the lowest lying states is not reproduced. In particular, 
the N(1440)Pii ('Roper') resonance and the first excited 
A, the P33(1600), which in the quark model belong to the 
N=2 oscillator shell, appear well below the states from the 
N=l shell. Furthermore, even the models with the fewest 
effective degrees of freedom predict many more states than 
have been observed, which is known as the 'missing reso- 
nance' problem. 

Since most states have been observed with elastic scat- 
tering of charged pious it is possible that the database 
is biased towards states that couple strongly to ttN. As 
an alternative, photon induced reactions, which nowadays 
can be experimentally investigated with comparable pre- 
cision as hadron induced reactions, have moved into focus. 
In order to avoid bias due to the resonance decay prop- 
erties, a large effort has been made during the last few 
years at tagged photon facilities to study many different 
final states [mil2] . So far, these experiments have mostly 
concentrated on the free proton. However, data from the 
neutron is also important, because it reveals the iso-spin 
composition of the electromagnetic excitation amplitudes. 
In extreme cases 'jNN* couplings may even be completely 
forbidden due to SU(6) selection rules [T?. Although due 
to the non-negligible spin-orbit mixing in the wave func- 
tions they are not strictly forbidden in more realistic mod- 
els, they remain suppressed and can be better studied us- 
ing neutron targets. Such experiments are of course com- 
plicated by the non-availability of free neutrons as targets, 
requiring coincident detection of recoil neutrons from light 
target nuclei and reaction models taking into account pos- 
sible nuclear effects on the observed cross sections. 

During the last few years the CBELSA/TAPS col- 
laboration has started in Bonn an extensive program for 
the study of quasi-free meson production reactions off the 



deuteron, including the nr] |14j . nr]' |15l . nut, nvr^Tr^, and 
nrjTT'^ final states. These are reactions with only neutrons 
and photons in the final state, which can only be inves- 
tigated with highly efficient electromagnetic calorimeters 
covering almost the full solid angle. In the present paper 
we discuss the experimental details and analysis proce- 
dures and summarize the results for the r/-photoproduction. 
Some results from this reaction have already been pre- 
sented in a letter [2], the results for the Ty'-channel have 
been published recently, and the results from the other 
channels will be published elsewhere. 



2 Quasi-free photoproduction of ?7-mesons 

Photoproduction of ry-mesons off the free proton has been 
previously studied in detail P^[r71[T51[ra|[^[?T|[^[^[^ 
[g51[ggl[g71[Ml[g g [gn i [gI] from the production threshold at 
ss707 MeV up to incident photon energies of 2.8 GeV. 
In the threshold region this reaction is completely domi- 
nated by the photoexcitation of the Sii(1535) resonance 
[22] . A detailed analysis of the angular distributions [TB] 
revealed a small contribution of the D13 (1520) via an inter- 
ference with the leading i?o+ -multipole of the Sn excita- 
tion. The effect is more pronounced for the beam asymme- 
try measured with linearly polarized photons |17ll26j and 
an analysis in the framework of the 'Eta-MAID' model 
[55] allowed the extraction of the tiny Nrj branching ra- 
tio (0.23±0.04 %) [Ml of the D13 resonance. In the range 
of the second Sn resonance, the Sii(1650), all models 
find a destructive interference between the two Sn states. 
The situation is less clear above this range, where dif- 
ferent analyses like Eta-MAID [33] and Bonn-Gatchina 
(BoGa) [3S] propose different resonance contributions (see 
e.g. [26]). 

So far, the isospin degree of freedom was almost ex- 
clusively explored in the excitation range of the first Sn 
resonance. The results of quasi-free production off the 
deuteron and ^He with and without detection of recoil 
nucleons [3M371I351I3U] are all consistent with: 



(1) 



where cr„ and ap are the cross sections for yy-photoproduction 
off neutrons and protons, respectively. The resulting ratio 
of the helicity couplings A"^,^ and A^,^ for the Sii(1535) 

is (n] 

K/2l/l^?/2l= 0-82 ±0.02 . (2) 



The investigation of coherent 77-photoproduction off the 
deuteron [57[|iD] . ^He [iT, and ''He [3^, and the compari- 
son of the interference terms in the angular distributions of 
quasi-free production off the proton and the neutron [35] 
have shown that the Sii(1535) excitation is dominantly 
iso- vector, so that [T2] 



^i/2/^1/2 =0.09 ±0.01 
where A^^% is the iso-scalar part of the amplitude. 



(3) 
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Fig. 1. Comparison of measured (Jn/ap ratios from |38II39| to 
model predictions. Eta-Maid model [3^ (MAID 1 (solid curve): 
full model, dotted curve: full model after folding with momen- 
tum distributions of bound nucleons, MAID 2 (dashed curve): 
only Sii(1535) and background). Dash-dotted curve: chirally 
coupled channel dynamics, Kaiser et al. [421I43J . E*: effective 
photon energy corrected for Fermi momentum. 



Previously measured and predicted cross section ra- 
tios are summarized in Fig. [T] The Eta- MAID model [55] 
agrees well with the 2/3 ratio in the Sii(1535) range, but 
predicts for higher incident photon energies a significant 
rise due to the contributions from other resonances. This 
effect should still be visible for quasi-free cross sections 
measured for nucleons bound in the deuteron, which is 
demonstrated by the dotted curve in Fig. [T] The largest 
contribution to the rise comes in this model from the 
Di5(1675) resonance, which is one of the states which due 
to the Moorehouse selection rules [l3j should have much 
larger electromagnetic couplings to the neutron than to 
the proton. However, the Nrj branching ratio of this state 
was determined in the model from a fit to the proton data 
as 17 %. This rather large value is in confiict with other 
results; PDG [31] gives only an upper limit of < 1%. 

A strong rise of the ratio to higher incident photon 
energies was also predicted in the framework of a chiral 
coupled channel model by Kaiser and collaborators [HI 
H5] . But in this case the rise is not related to the con- 
tribution of higher lying resonances, it is related to the 
properties of the Sii(1535), which in this model is not 
interpreted as a genuine qqq state but as a dynamically 
generated K U-molecule like state. 

Finally, also in the framework of the chiral soliton 
model [HIHS] a state was predicted in this energy range, 
which has a much stronger photon coupling to the neutron 
than to the proton and a large decay branching ratio into 
Nrj. This state is the nucleon-like member of the conjec- 
tured anti-decuplet of pentaquarks, which would be a Pn 
state. Exact SU(3)i? would forbid the photo-excitation of 
the proton to the proton- like member of the anti-decuplet. 
But even after accounting for SU(3)f violation the chiral 



soliton model predicts 44 that the photo-excitation of 
this state is suppressed on the proton and should mainly 
occur on the neutron. Kim et al. [JB] have calculated the 
magnetic transitions moments for the anti-decuplet states 
and found a considerable enhancement for the excitation 
of the nucleon-like state on the neutron with respect to 
the proton. 

Strong efforts have recently been undertaken at various 
facilities (GRAAL in Grenoble [IZj, ELSA in Bonn [H] 
and LNS in Sendai [4S1) to extract reliable results for the 
jn — >■ nr] reaction at higher incident photon energies. The 
somewhat unexpected finding in all experiments is not 
only a significant rise of the cross section on the neutron 
with rising incident photon energy, but a fairly narrow 
peak in it, which has no counterpart in the reaction off 
the proton. 

The nature of this structure is still unknown and many 
different suggestions have been put forward in the litera- 
ture. They include various types of coupled channel ef- 
fects involving known nucleon resonances [491150] or the 
opening of production thresholds 51 , but also scenarios 
with contributions from intrinsically narrow states have 
been discussed [5^155] . Very recently, Kuznetsov et al. [SI] 
have reported a narrow structure in Compton scattering 
off the quasi-free neutron with a similar mass and width, 
although not with large statistical significance (w4.6(t). 
Such an observation would make explanations with com- 
plicated interference and threshold effects less likely since 
there is no good reason why they should appear similar for 
such different reactions. A nucleon resonance with strong 
photon coupling to the neutron, on the other hand, nat- 
urally would be expected to show up also in Compton 
scattering off the neutron. 

In this paper, we will summarize and discuss in detail 
the results from the ELSA experiment reported partly in 
|14) and in addition present a new analysis which removes 
the effects of Fermi motion from the data through a kine- 
matic reconstruction of the involved nucleon momenta. In 
this way, more stringent constraints can be put on the 
intrinsic width of the peak-like structure. 

The paper is organized in the following way. The exper- 
imental setup is described in Sec. [3] Details of the analysis 
are discussed in Sec. |4l This includes the calibration of all 
detector components, the identification of photons, recon- 
struction of mesons, and detection of recoil nucleons, the 
absolute normalization of the cross sections and a thor- 
ough discussion of systematic uncertainties. The results 
are summarized in Sec. [5] We first compare the data to 
previous results from quasi-free photoproduction off the 
deuteron and the results for the quasi-free proton to free 
proton data. In the following two subsections the results 
for quasi-free photoproduction are first discussed in depen- 
dence of the incident photon energy. In this analysis the 
width of narrow structures is dominated by nuclear Fermi 
motion. In a second analysis results are constructed in de- 
pendence of the nucleon - meson invariant mass in the final 
state calculated from the four vectors of the observed par- 
ticles. In this way, Fermi motion effects are removed and 
only instrumental resolution must be considered. 
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Fig. 2. Overview of the experimental setup. The electron beam enters from the left side. Scattered electrons were detected in the 
counters of the tagging spectrometer (cf. Fig. (3]). The liquid deuterium target was mounted in the center of the Crystal Barrel. 
The forward range was covered by the TAPS detector. The Time-of- flight wall was mounted, but not used in the experiment. 
Beam intensity was monitored at the end of the beam line. 



3 Experimental setup 

The experiment was performed at the electron stretcher 
accelerator facility ELS A in Bonn [551156] . which can de- 
liver electron beams with intensities of a few nA for en- 
ergies up to 3.5 GeV. The overall setup of the experi- 
ment is shown in Fig. |5J The incoming electron beam is 
impinging on the radiator mounted in a goniometer. For 
the measurement discussed here, electron beam energies 
of 2.6 GeV and 3.2 GeV were used. The largest part of the 
3.2 GeV beam time was done with a copper radiator foil of 
0.3 % radiation lengths thickness, producing unpolarized 
bremsstrahlung. For the rest of this beam time and for the 
beam times with 2.6 GeV electron energy a diamond ra- 
diator was used to produce coherent bremsstrahlung with 
maximum linear polarization around 1 GeV incident pho- 
ton energy (see [571 for details about linear polarization 
of the photon beam) . 

Electrons which have emitted bremsstrahlung, are then 
deflected downwards by a dipole magnetic field onto the 
focal plane of the tagging system, where energy and tim- 
ing information is extracted. Unscattered electrons are 
stopped in the beam dump. The bremsstrahlung photons 
are almost co-linear with the incident electron beam, pass 
through a hole in the magnet yoke and irradiate the liq- 
uid deuterium target mounted in the center of the Crystal 
Barrel detector. The core of the detector system of the 
tagging spectrometer (see Fig. [3]) are 14 partially overlap- 
ping plastic scintillator bars each with photo-multipliers 
at both ends. This system covers 22 % to 92 % of the in- 
coming electron beam energy Eg . A second layer has been 
added in front of the bars for better energy resolution. 
A scintillating fiber detector (2 x 240 fibers arranged in 
two layers to partly overlap) covers photon energies from 
18 % to 80 % of Eo- This defined the maximum tagged 
photon energy of 2.5 GeV. The fiber detector provides an 
energy bin width of «1.5 % for the lowest incident pho- 
ton energies and «0.1 % at the high energy. In principle, 
the range from 80 % to 92 % of Eo can be tagged with 
an additional wire chamber. The chamber was, however, 
not used in the present experiment since due to the small 
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Fig. 3. Setup of the tagging spectrometer. Unscattered elec- 
trons stop in the beam-dump. Scattered electrons pass a two 
layer detection system. Scintillating fiber detectors (good po- 
sition resolution) and scintillator bars (good time resolution) 
cover photon energies up to 80 % of the electron beam energy. 
The part covered by a multiple wire chamber (MWPC) (low 
energy electrons corresponding to photon energies above 80 % 
of the electron beam energy) was not used in the experiment. 



cross section the statistical quality of data in this range 
was not sufficient. 

The target consisted of a vertically mounted cryostat 
attached to a tube entering the Crystal Barrel detector 
from the upstream side. The target cell itself was a cap- 
ton cylinder (0.125 mm foil thickness) with a diameter of 
3 cm and length of 5.275 cm, filled with liquid deuterium 
(surface thickness 0.26 nuclei/barn). 

Reaction products emerging from the target have been 
detected with a combined setup (see Fig.|3|) of the Crystal 
Barrel detector [55] and the TAPS detector [5^1150] . In the 
configuration used, the Crystal Barrel consisted of 1290 
Csl (Tl) crystals of 16 radiation lengths Xg all mounted 
in a target pointing geometry. It covered the full azimuthal 
angle for polar angles between 30° and 168°. The forward 
angular range was covered by the TAPS detector [5^160] . 
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Fig. 4. Arrangement of the Crystal Barrel and TAPS detec- 
tors. Upper part: side view, lower part: front view of the TAPS 
wall: left hand side: logical segmentation for the LED-low trig- 
ger, right hand side: logical segmentation for the LED-high 
trigger (see text). 



This component was made of 528 BaF2 crystals of hexag- 
onal shape with an inner diameter of 5.9 cm and a length 
of 25 cm corresponding to 12 Xq. They were arranged in 
a wall-like structure as shown in the lower part of Fig. 21 
covering polar angles down to 4.5°. The front face of the 
BaF2 wall was located 1.18 m from the center of the tar- 
get. 

The two calorimeters have a comparable energy reso- 
lution [53160] 

a^_, 2-3% ^^^ 



E 



^E/GeY 



Since the impact points of photons are determined from 
the center of gravity of the electromagnetic showers, the 
angular resolution is better than the granularity of the 
crystals. It is 1.5° (tr) for the CB (52 for photons with 
energies above 50 MeV and 1.25° in TAPS. Angular reso- 
lution for recoil nucleons, which do not produce extended 
clusters, is closer to the granularity. 

Both parts of the calorimeter were equipped with de- 
tectors for the identification of charged particles. A three- 
layer scintillating fiber detector (Tnner'-detector) [ST] was 
mounted inside the Crystal Barrel around the target, cov- 
ering polar angles between 28° and 172°. The outer layer 
(diameter 12.8 cm, 191 fibers) runs parallel to the z-axis. 
The middle layer (0= 12.2 cm, 165 fibers) is wound anti- 
clockwise at an angle of 25° with respect to the z-axis and 



the inner layer (0= 11.6 cm, 157 fibers) is wound clock- 
wise at 25°. All fibers have 2 mm diameter. This orien- 
tation allows the reconstruction of the spatial coordinates 
of the intersection point of the charged particle trajec- 
tory with the detector (see [55] for details). The TAPS 
detector was complemented with a charged-particle-veto 
(CPV) detector built of 5 mm thick plastic scintillators 
of hexagonal shape and the same dimensions as the front- 
face of the BaF2 crystals, so that each detector module 
had its individual veto detector. The CPV was read out 
by wave-length shifting fibers connected to multi-anode 
photo- multipliers. 

The fast BaF2 modules of the TAPS detector were 
read out by photo-multipliers, but the Csl crystals of the 
Barrel were read by photo-diodes, which do not provide 
any time information. This had important consequences 
for the hardware trigger. Signals for the first-level trig- 
ger could only be derived from the TAPS detector. For 
this purpose each module was equipped with two indepen- 
dent leading-edge discriminators, combined in two differ- 
ent ways into logical groups (see Fig.S]). For most of them 
(ring 12-5 from outer edge to center) a lower threshold 
was set to w55 MeV (LED-low). It was set to 80 MeV, 
135 MeV, 270 MeV for rings 4, 3, 2, respectively. The 
inner-most ring was not used in the trigger. The LED-high 
thresholds were set to 70 MeV for rings 9-7, rising from 
105 MeV (ring 6) to 180 MeV (ring 2). Again, the inner- 
most ring was not used in the trigger and the three outer 
rings (block G) had no leading edge discriminators for the 
high threshold. The first-level trigger then included two 
components: at least two LED-low discriminators from 
different logical sectors above threshold, or at least one 
LED-high discriminator above threshold. In the second 
case, a second-level trigger from the FAst Cluster Encoder 
(FACE) of the Crystal Barrel, indicating at least two sep- 
arated hits in the Barrel, was required in addition. All first 
level triggers thus required detection of one or two pho- 
tons in TAPS, which covered only a small part of the solid 
angle. Therefore, only reactions with relatively high pho- 
ton multiplicity like ^d — >■ npn'^Tr^ — >■ npA'y, jd — > npr] — > 
npSn'^ — )> np6^, or jd -^ npr]' — ^ nprj2Tr^ — ?> np6j could be 
recorded. Two-photon decays of single meson production 
reactions like ^d -^ npir'^ — >■ np2j or jd — >■ npi] — >■ np2^ 
were not taken. It should be mentioned that this restric- 
tion of the detector setup applies only to reactions off the 
(quasi)-free neutron. In measurements off the free pro- 
ton, the recoil protons provide trigger signals from the 
Inner-detector (or TAPS). The TAPS modules were ad- 
ditionally equipped with constant fraction discriminators 
(CFD), with thresholds around 10 MeV, which were used 
for the high resolution measurement of the time and gen- 
erated the read-out pattern of the detector matrix. 

The last component of the setup is the so-called 7- veto 
detector at the end of the photon-beam line, which con- 
sisted of nine lead-glass crystals, and counted the photons 
which passed through the target without interaction. It 
was used to monitor the time dependence of the photon- 
beam intensity. Furthermore, it was used to measure the 
tagging efficiency, i.e. the fraction of photons, correlated 
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with electrons, which impinged on the target. This was 
done in special tagging efficiency runs, where at reduced 
beam intensity the tagger was used as trigger and the 
number of scattered electrons counted in the tagger was 
compared to the number of coincident photons impinging 
on the 7- veto. 



4 Data analysis 

In this chapter we discuss the different analysis steps of the 
experiment, from the most basic calibration procedures 
for the different detector components, over the identifi- 
cation of photons and particles, the identification of spe- 
cific reaction channels, to the absolute normalization of 
cross sections. Although in this paper only final states 
with ry mesons are discussed, a large part of this chapter 
applies also to the other, simultaneously measured reac- 
tions, which will be reported elsewhere. Therefore we dis- 
cuss in this paper the analysis steps in some detail also in 
view of possible systematic uncertainties. A full account 
of all analysis procedures is given in [^5] . 



4.1 Calibration procedures 

4.1.1 Tagging system 

The tagging system has two tasks: the event- by- event def- 
inition of the incident photon energy via a coincidence be- 
tween the focal plane detectors and the reaction detector 
and the monitoring of the photon flux from the counting 
of the deflected electrons. 

For the present analysis energy and timing information 
was obtained with the scintillating fiber detector. Since for 
any scattered electron the bending radius p in the mag- 
netic field B is related to its momentum p by 



P = 



P 
Bq 



(5) 



where q is the electron charge, its energy follows from the 
position in the focal plane detector, i.e. the number of 
the responding fiber. The energy calibration of the scin- 
tillating fiber detector was done in two steps. As a start- 
ing point a polynomial calibration function was calculated 
from the measured field map of the dipole magnet and the 
positions of the scintillating fibers. Initially, this calibra- 
tion was done for an electron beam energy of 3.2 GeV. 
For other beam energies, the magnetic field was adjusted 
such that the unscattered electrons always followed the 
same trajectory so that the calibration function could be 
simply scaled by the electron beam energy. This calibra- 
tion was checked by direct injection of a very low intensity 
electron beam with removed radiator. At a fixed magnetic 
field setting for 3.2 GeV (5=1.413 T) ELSA beams with 
four different energies (1.2, 1.5, 2.0, 2.5 GeV) were used. 
A further cross check was done with the position of 
the peaks from coherent bremsstrahlung produced in a 
diamond crystal, which was used for measurements with 
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Fig. 5. Relative TAPS-tagger time spectrum for all events 
with photons in TAPS (solid histograms) and for events with 
r;-decay photons in TAPS (dashed histograms). The 2 ns time 
structure of the beam is reflected in the fine structure of the 
accidentals. The shaded areas indicate the prompt (A) and 
random (B, C, D) coincidence areas (see text). 



linearly polarized photon beams. The systematic uncer- 
tainty of this calibration is on the order of the energy bin 
widths of the fiber detector. 

The measurement of the relative timing between the 
scattered electrons and the reaction products was done 
with the focal plane scintillating fiber detector and the 
TAPS detector. The modules of both detector systems are 
equipped with individual TDC's. In case of the focal plane 
detector the start signal came from the individual fibers 
and the stop from the trigger signal. The TAPS TDC's 
were started with the trigger signal and stopped by the 
individual CFD signals. The time calibration of the tagger 
TDC's was 64 ps/channel. After alignment of all channels 
a time resolution of 1.6 ns (FWHM) was achieved for the 
sum spectrum, which is shown in Fig. [S] The average mul- 
tiplicity of hits in the focal plane detector per event in the 
coincidence time window of 460 ns was ~20, resulting in a 
significant random background. After a cut on the prompt 
peak (±3 ns), the remaining background corresponding 
to region B in Fig. [5] was removed from all shown spec- 
tra in the usual way by subtraction of the events from 
the areas D and C normalized by the ratio of the areas 
B/(D-|-C). The random background was much less impor- 
tant for events with photons from an identified 77-meson 
(dashed histograms in Fig. [5]). Also the small asymmet- 
ric tail of the prompt peak (cf. Fig. [5]), due to particles 
misidentified as photons, vanishes in this case. 
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4.1.2 TAPS 

A precise calibration of the time measurement with TAPS 
was not only important for the suppression of random 
background. It was also the basis of the measurement of 
the kinetic energy of recoil nucleons with the time-of- flight 
method. The slightly varying gaines of the TDC's were 
measured for each channel with pulser signals of known 
delay, fed into the electronic chain. The alignment of all 
channels was done with events with two or more photons 
detected in TAPS. For this procedure, photons were iden- 
tified with the help of the veto detector and an invari- 
ant mass analysis, which accepted only decay photons of 
7r° mesons. The resulting TAPS - TAPS coincidence time 
spectrum had a resolution of 650 ps (much better than the 
TAPS - tagger timing) and a peak-to-background ratio of 
100. 

The energy calibration was done in three steps. In the 
first step the energy deposition of minimum ionizing cos- 
mic myons («38.5 MeV) were used for a relative calibra- 
tion of all detector modules. Such a calibration does not 
take into account differences in the detector response to 
the energy deposition of myons and electromagnetic show- 
ers as well as shower leakage. Therefore, in the second step 
the linear term of the calibration was adjusted for all mod- 
ules in an iterative procedure with the invariant mass peak 
of TT^ mesons. In the third step, a quadratic correction of 
the energy calibration was introduced using the position 
of the invariant mass peak of the rj meson. A final accu- 
racy of ±1 % for the position of the 7r°- and ry-peaks was 
achieved. More details can be found in 1631. 



4.1.3 Crystal Barrel 

Due to the read-out with photo-diodes, the CB delivers 
only energy information. A calibration with cosmic myons 
is not practical, due to the different geometries and orien- 
tations of the individual modules. The digitization of the 
detector signals is done for two different dynamic ranges. 
The calibration for the low-gain chain was again done with 
an iterative procedure using the tt'^ invariant mass peak. 
The high gain branch was calibrated by the injection of 
laser light of known intensity into the crystals. 



4.2 Identification of particles and reaction channels 

Electromagnetic showers, depending on their energy, will 
in general produce signals in an extended 'cluster' of scin- 
tillator modules. In a first step of the analysis all hits in 
the calorimeters were grouped into 'clusters' of adjacent 
crystals and the energy sums and geometrical centers of 
gravity of the 'clusters' were extracted (see [52] for de- 
tails). In the next step, the TAPS- veto detector and the 
Inner-detector were used to separate neutral hits (photon 
and neutron candidates) from charged hits (proton can- 
didates) in the detectors. In TAPS a hit was assigned to 
'charged' if the veto of any module from the cluster or the 
veto of any neighbor of the central module of the cluster 



had fired (even if the neighbor module itself had no sig- 
nal above threshold). The latter condition is important for 
charged particles with relatively large impact angles which 
may traverse the edge of a veto but deposit their energy 
in the neighbor module. All other hits were assigned to 
'neutral'. In the Barrel, a hit was assigned to 'charged', 
if at least two layers of the Inner-detector had recorded 
a hit within ±10° of azimuthal angle. The efficiency of 
the Inncr-dctcctor for this condition is 98.4 % [5T]. It was 
assigned to 'neutral' if no layer had fired within this az- 
imuthal angle, which results in a probability of « 0.04 % 
to misidcntify a charged particle as neutral hit. Events 
with hits with one responding layer of the Inner-detector 
were discarded. 

After these assignments three partly overlapping classes 
of events were included into the analysis. Class (a) in- 
cluded events with exactly six neutral hits and exactly one 
charged hit. Class (b) included events with exactly seven 
neutral hits and no charged hit (in this case was addition- 
ally required that the Inner detector had not fired at all) . 
Class (c) included all events from classes (a) and (b) and 
in addition all events with exactly six neutral hits and no 
charged hit. Class (a) corresponds to quasi-free production 
off the proton jd -^ {n)pr] with coincident detection of the 
recoil proton. Class (b) corresponds to quasi-free produc- 
tion off the neutron jd -^ {p)nri with coincident detection 
of the recoil neutron. Class (c) corresponds to the inclu- 
sive reaction ^d -^► {np)rj, without any condition for the 
recoil nucleon (may be detected but is not required) . 

These events were subjected to a combined invariant 
and missing mass analysis. In the first step of the invari- 
ant mass analysis the invariant masses of all combinations 
of three disjunct pairs of neutral hits were calculated. In 
the case of six neutral hits (events with proton candidate 
or without candidate for recoil nucleon) these are 15 dif- 
ferent combinations among which the 'best' combination 
was chosen by a x^-test, minimizing 



771 Tfo)^ 



^ (Zimfe (77))2 



(6) 



for all disjunct combinations where m^o is the pion mass 



and the 






are the invariant masses of the photon pairs 



with their uncertainties Z\m^^ calculated for each pho- 
ton pair from the energy and angular resolution of the 
detector. For events with seven neutral hits (events with 
neutron candidate) one must in addition loop over the un- 
paired hit giving rise to 105 combinations. Once the 'best' 
combination was determined, in all cases a cut between 
110 MeV - 160 MeV was applied to the invariant masses. 
Only events where all three pairs of the best combination 
passed this cut were kept. Subsequently, for events with 
seven neutral hits, the residual hit was taken as neutron 
candidate. 

The nominal mass of the pion was then used to improve 
the experimental resolution. Since the angular resolution 
of the detector is much better than the energy resolution 
this was done by re-calculating the photon energies and 
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Fig. 6. Invariant (left hand column) and missing mass (right 
hand column) spectra for 7;-mesons in coincidence with recoil 
neutrons for three ranges of incident photon energy. Invariant 
masses: shaded (blue) signal after missing mass cut (scaled up 
by indicated factors). Dashed lines: applied invariant mass cut. 
Missing mass data: open symbols represents data for indicated 
cut on invariant mass. Black dots: background subtracted by 
fitting invariant mass spectra for each bin of missing mass. 
Simulations: dashed (dotted, dash-dotted) curves: simulation 
of ?7 (rin, rj') final states. Solid (red) curves: sum of simulations. 
Shaded (blue) areas: accepted events. 



momenta from the approximation: 



^1,2 — -^1:2 



(7) 



where Ei^2 are the measm'ed photon energies, £[ 2 the 
re-calculated energies, m^o is the nominal tt^ mass, and 
rrijj the measured invariant mass. The slow E~^'^ energy 
dependence of as/E from Eq. Uhas been neglected here. 

In the last step, the invariant mass of the six hits as- 
signed as photons was calculated. The result is shown in 
Fig. [S] (left hand side) for three different ranges of incident 
photon energy. The figure shows the most difficult case 
for events with neutron candidates, where combinatorial 
background in the x^ analysis of the best combination is 
higher than for events with six neutrals only. 

After the invariant mass analysis all events were sub- 
jected to a missing mass analysis, where the recoil nucleon, 
no matter if a candidate was detected or not, was treated 
as missing particle. The missing mass Am of the reaction 
was calculated for quasi-free production of 77 mesons off 



Am 



N 



mN 



(8) 



Pj, are the four- 



where ttin is the nucleon mass, P^, Pn, 
momenta of the incident photon, the initial state nucleon 
(assumed to be at rest), and the produced 77- meson. The 
resulting distributions peak around zero. They are some- 
what broadened by the momentum distribution of the 
bound nucleons, which was neglected. Typical results are 
shown on the right hand side of Fig. [S] The open symbols 
correspond to the data after the cuts on invariant mass 
shown on the left hand side of the figure have been ap- 
plied. At higher incident photon energies background in 
particular from rjiT final states is visible, where the pion 
has escaped detection. A smaller background at the high- 
est incident photon energies is due to the 77' -^ tjit'^tt'^ de- 
cay. Cleaner results are obtained when the invariant mass 
spectra are generated for each bin of missing mass and 
fitted with the 77 line shape and a background polyno- 
mial. This analysis demonstrates (compare open and filled 
symbols in Fig. [HI right hand side) that the background 
underneath the 77 invariant mass peak contributes mainly 
to large missing masses. Only a small contribution from 
triple tt" production appears in the missing mass peak re- 
gion. The shape of the missing mass structures for the 77, 
777r, and rj' final states has been generated with a Monte 
Carlo simulation (see below) . The sum of these contribu- 
tions reproduces the measured spectra. 
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Fig. 7. Simulated distribution of recoil nucleons depending 
on incident photon energy and rj cm polar angle. Nucleons 
corresponding to the (blue) shaded area are emitted into the 
solid angle covered by TAPS, nucleons from the white area into 
the Barrel. Due to Fermi smearing the borders are not sharp. 



The final analysis was done in the following way. For 
incident photon energies below 800 MeV, where the reac- 
tion jd -^ tj-kX is kinematically forbidden and the miss- 
ing mass spectra are background free, the signals were 
obtained by fitting the invariant mass spectra with peak 
shape and background polynomial. At higher incident pho- 
ton energies, where background is no longer negligible. 
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first a cut on missing mass was applied, which only ac- 
cepted events on the left hand side of the peak (blue 
shaded areas in Fig. [51 right hand side). The corresponding 
invariant mass spectra are shown on the left hand side of 
the figure as (blue) shaded histograms. The small residual 
background underneath these invariant mass peaks was 
fitted and removed. This analysis was done for each bin of 
incident photon energy and rj center-of-momentum (cm) 
polar angle (0*). The missing mass cut was chosen very 
restrictive, at the price of the loss of half the counting 
statistics. This was done to exclude any background con- 
tamination, which could create artificial structures in the 
excitation function around 1 GeV where the background 
in the missing mass spectra starts to appear. 

Finally, we discuss the assignment of recoil nucleons 
to the events. Fig. [7] shows which recoil nucleons need to 
be detected in TAPS and which in the Barrel, depend- 
ing on the reaction kinematics. For incident photon en- 
ergies below « 850 MeV all recoil nucleons are emitted 
into the solid angle covered by TAPS. At higher energies 
only nucleons corresponding to backward emission of the 
ry-mesons (approximately cos(6'*)<-0.2) are detected in 
TAPS. 



removed by the missing mass cut since they originate from 
r]TT final states. 

In case of the TAPS detector, which was placed 1.18 m 
downstream from the target, an independent check of the 
recoil nucleon identification was done via a time-of-flight 
versus energy analysis of the recoil particles. This is illus- 
trated in Fig. [H In the spectrum of charged recoil particles 
detected in coincidence with an 77-meson identified by in- 
variant mass (upper left corner) a clear proton band is 
visible. After the missing mass cut (lower left corner) al- 
most all background is removed. The spectrum of neutron 
candidates after application of 77 invariant and missing 
mass cut (lower right corner) does not show any significant 
trace of the proton band. The maximum contamination of 
the TAPS neutron sample by recoil protons was estimated 
from these spectra to be below 3 % (Monte Carlo simula- 
tions indicate a maximum contamination at the 4% level) . 
Finally, for the rin" final state (upper right corner), which 
is not further discussed in the present paper, also a clear 
band for recoil deuterons is visible. This band is missing 
for the 77 final state since coherent production of single rj 
mesons off the deuteron is strongly suppressed [ID] . 




t[ns] 



Fig. 8. Time-of-flight versus energy spectra measured with 
TAPS for different event types. Upper left: 77 — >■ 67 plus one 
charged particle, upper right: 77 — > 27 and n° — >■ 27 plus one 
charged particle, lower left: r; — > 67 plus one charged particle 
Cp') and missing mass cut, lower right: 77 — >■ 67 plus one neutral 
particle {'n') and missing mass cut. 



Nucleons detected in the Barrel are identified as pro- 
tons or neutrons as discussed above with the help of the 
Inner-detector. Due to the three-layer structure of this de- 
vice the probabilities to misidentify a proton as neutron 
or vice versa are negligible. Charged pions might be also 
misidentified as protons. However, such events are reliably 



4.3 Determination of cross sections and systematic 
uncertainties 

The extraction of cross sections from the analyzed yields 
requires several pieces of information: the 77 -^ 67 decay 
branching ratio, the target surface density, the incident 
energy dependent photon flux, and the detection efficiency 
of the combined Crystal Barrel - TAPS calorimeter with 
particle identification detectors for photons, protons, and 
neutrons including all analysis cuts. 

The decay branching ratios for 77 — > Stt'' and tt^ -^ 2j 
are taken from the Particle Data Group [33] as (32.56± 
0.23)% and (98.823± 0.034)%, respectively, resulting in 
a total branching ratio of 31.42 % with a negligible sys- 
tematic uncertainty. The target surface density (5.3 cm 
long liquid deuterium target, density p ~ 0.169 g/cm^) 
was 0.26 deuterons/barn with a systematic uncertainty of 
«2%. 

The determination of photon flux and detection effi- 
ciency are discussed in detail in the following sub-sections. 



4.3.1 Flux normalization 

The data have been taken in five blocks of beam time 
with different electron beam energies and settings of lin- 
ear beam polarization, which we label A, B, C, D, E. They 
are summarized in Tab.[T] Due to these settings the statis- 
tical quality of the data is less good above 2 GeV incident 
photon energy. 

The photon flux was determined in the following way. 
The number of scattered electrons detected by the tagger 
focal plane counters (scintillating fibers) was determined 
with scalers. The scalers were not life time gated, but the 
experiment dead time (typically 40 %) was recorded with 
separate scalers and corrected. Due to an electronics/data 
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Table 1. Summary of beam times. E^- : electron beam energy, 
E^^: maximum energy of tagged photons, Epd'. energy of max- 
imum linear photon beam polarization, (Po'- energy integrated 
electron flux on the tagging detectors. Total life time: beam 
time multiplied by acquisition life time fraction. 



characteristics 


A 


B 


C 


D 


E 


E,- [GeV] 


2.6 


2.6 


3.2 


3.2 


3.2 


E^, [GeV] 


2.0 


2.0 


2.5 


2.5 


2.5 


Epoi [GeV] 


1.0 


1.0 


unpol. 


1.2 


1.6 


total life time [h] 


138 


18 


189 


25 


25 


'Po [10'e-/s] 


1.75 


1.6 


1.6 


2.8 


2.8 



acquisition problem the scaler information was incorrectly 
handled for part of the data. This resulted in a loss of the 
absolute normalization for the beam time blocks B, D, and 
E, which, however, represent only a minor part of the total 
statistics. In order to recover the total available statistics, 
these three data blocks were relatively normalized to the 
data sets A, C, using the total cross sections for inclusive 
7y-production and double 7r° production. 

The absolute normalization furthermore requires the 
measurement of the fraction of correlated photons that 
impinge on the target. This was done with special tagging 
efficiency runs, where the intensity of the photon beam 
was directly measured at the end of the beam line at re- 
duced electron beam intensity. The average tagging effi- 
ciency, which is rather flat as function of beam energy, was 
«95 % for beam time A and «75 % for beam time C (for 
the latter the beam quality was less good). 
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Fig. 9. Main figure: incident photon flux for beam time A 
(red dots) and beam time C (open black circles). Insert: ratio 
of inclusive 77-photoproduction cross sections for the two beam 
times. 



The photon flux, i.e. the product of the spectrum of 
electron scalers throughout the tagger detector and the 
tagging efficiency, for beam times A and C are compared 
in Fig. ini The overall shapes follow the typical 1/E~^ be- 



havior of bremsstrahlung. For beam time A the peak from 
coherent bremsstrahlung is visible around I GeV photon 
energy. The visible structures are due to systematic ef- 
fects of individual tagger counters. They are also present 
in the extracted yields and cancel in the cross sections. 
The insert shows the ratio of the total cross sections for 
inclusive 7y-photoproduction obtained with the respective 
photon fluxes for beam times A and C. For most of the 
energy range systematic deviations are between % and 
-10 %. At the lowest incident photon energies they are 
around -1-10%. For the measurement with the 3.2 GeV 
electron beam, this energy region corresponds to the very 
edge of the tagged range, where systematic uncertainties 
tend to be larger. For the final result all data sets were 
averaged according to their statistical weights and a sys- 
tematic photon flux uncertainty of 10 % was estimated. 



4.3.2 Monte Carlo simulations and detection efficiency 

The detection efficiency of the Crystal Barrel/TAPS setup 
was modelled with Monte Carlo simulations based on the 
GEANT3 package [M]. The simulations include all rel- 
evant properties of the calorimeter, including geometri- 
cal acceptance, charged particle identification, trigger ef- 
ficiency, response of all detector modules, and analysis 
cuts. They include also information about inefficient or 
malfunctioning individual detector modules. This means 
that the extracted detection efficiencies are effective ones, 
which cannot be directly applied to data sets taken under 
different conditions. As far as photon detection and iden- 
tification of 7/-mesons via the ry — > Stt'^ — > 67 decay chain 
are concerned, some details are already given in |66j . which 
used an identical setup for the study of 77-photoproduction 
from heavy nuclei. However, here in addition the detection 
efficiency for recoil nucleons plays a crucial role. The sim- 
ulations were done with the GEANT-CALOR program 
package [67j . which is optimized for hadronic interactions 
from the few MeV to several GeV range, including the 
interactions of low energy neutrons. 

The results for the recoil nucleon detection efficiencies 
are summarized in Fig. 1101 For both, neutrons and pro- 
tons, in the first step isotropically distributed nucleons 
were simulated. The detection efficiency was determined 
as function of kinetic energy T and laboratory polar an- 
gle 0. As expected, angular dependencies of the efficiency 
(not shown) are negligible for the CB apart from the very 
edge of the detector. They are more important in the 
TAPS detector (changing angle of incidence towards the 
outer edge, increasing discriminator thresholds towards 
the beam pipe). Fig. [TU] (black open circles) shows the 
kinetic energy dependence of the detection efficiency ob- 
tained with this simulation, averaged over the polar angle 
range of the detectors. 

In case of the proton, the simulation can be checked 
with data from the liquid hydrogen target measured with 
the same setup. For this purpose, the reactions "fp -^ 
pr] -^ p6j and jp — ?• p7r°7r° -^ pA'y were analyzed using 
similar procedures for invariant and missing mass anal- 
yses as discussed above. The proton detection efficiency 
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Fig. 10. Detection efficiency for recoil nucleons. Upper row: 
neutrons, lower row protons, left hand side: Crystal Barrel, 
right hand side: TAPS. For all figures: ei (open black cir- 
cles): MC simulation for isotropically emitted nucleons with- 
out further detector hits, £2 (blue squares): simulation with 
six decay photons from r;-meson decays, ea (red dots): exper- 
imentally determined proton detection efficiency from rj and 
27r° photoproduction from proton (hydrogen) target (see text). 
Solid (red)line: measured neutron detection efficiency for CB 
at CERN [BS]. 



ep(T, 0) is then simply given as the ratio of events with 
detected proton to all events of the respective reaction. 
The result, again averaged over the polar angles, agrees 
quite well with the simulation (see Fig. [TUl red dots). In 
case of the neutron, experimental information comes from 
the measurement of the neutron detection efficiency of the 
CB when it was installed at LEAR at CERN [55]. There- 
fore, in this case the simulation was done with the pa- 
rameters (e.g. thresholds) characterizing the CB setup at 
CERN. The results are also compared in Fig. [TU] Agree- 
ment is quite good at larger T. Some discrepancies are 
visible at low T, however, here the result is strongly de- 
pendent on detector thresholds and neutron energy cal- 
ibration, which may not have been reproduced perfectly 
in the present simulation of the old CB setup. In case of 
TAPS the neutron detection efficiency had been measured 
in an experiment at the Mainz MAMI accelerator using 
the 7p -^ riTT^TT"'" reaction for kinetic energies below 250 
MeV [38]. The results are consistent with the simulation 
when the conditions of the Mainz setup are used (at 250 
MeV simulated: ei « 18.5%, from data: ei « 19.1%). 

Finally, it must be considered that neutrons from the 
7n — >■ nrj reaction are identified in this experiment in 
events with seven neutral hits by first assigning six hits 
via the invariant mass analysis to the 77 — > Stt*^ -^ 67 
decay chain. This introduces efficiency losses due to com- 
binatorial background. Therefore, a simulation was done, 
where events from 771 — )> nq were produced with a phase 
space event generator. The analysis then mimicked the 
whole reconstruction process, including the identification 
of the neutron out of the seven neutral hits. The result- 
ing 'effective' neutron detection efficiency is also shown in 



Fig. [TU] (blue squares). It is significantly lower than the 
'raw' neutron efficiency. The same kind of analysis was 
also done for the proton case. However, due to the addi- 
tional information from Inner- and TAPS-veto detectors, 
in this case the effects are small. As final result of this 
analysis, recoil detection efficiencies for protons and neu- 
trons are available as function of laboratory polar angle 
and kinetic energy. 

The detection efficiency for ry-mesons has been simu- 
lated as function of their laboratory polar angles and ki- 
netic energies in the same way as discussed in detail in |66) . 
The result is shown in Fig. [11] (left hand side) and com- 
pared to the distribution of the experimentally detected 
77-mesons (Fig. [11] right hand side). The branching ratio 
of the ?7 ^^ 67 decay is not included in this efficiency (i.e. 
shown is the efficiency to detect 77-mesons which decay 
with 100 % branching ratio into six photons). The ac- 
ceptance of the detector covers the full phase space of the 
reaction. However, the absolute values of the efficiency are 
small in particular for mesons at backward angles. This is 
due to the efficiency of the trigger, which was only based 
on photon detection in TAPS. Note that the recoil nu- 
cleon detection efficiencies discussed above do not include 
trigger efficiencies, since the meson is assumed to trigger. 




000^ 



G_[Oj 



Fig. 11. Left hand side: simulated detection efficiency (includ- 
ing trigger efficiency) for 77 — >■ 67 as function of r] laboratory 
kinetic energy and polar angle. Right hand side: distribution 
of observed rj mesons in dependence on the same parameters. 



The total detection efficiencies for the inclusive chan- 
nel (no recoil nucleons required), the reactions in coinci- 
dence with recoil protons, and in coincidence with recoil 
neutrons were then calculated in two different ways to es- 
timate systematic effects. 

In the first approach, 77-photoproduction off quasi-free 
nucleons was simulated with a phase-space event genera- 
tor, taking into account the effects of nuclear Fermi smear- 
ing. The simulated events were analyzed in the same way 
as the data and the detection efficiency tph was calculated 
from the ratio of simulated and detected events as func- 
tion of incident photon energy and meson cm polar angle. 
The results are shown for some ranges of incident photon 
energy in Fig. [1^ as open squares. This simulation is in so 
far model dependent as the correlation between meson en- 
ergies and angles as well as the correlation between mesons 
and recoil nucleons relies on the phase-space assumption. 
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Fig. 12. Total detection efficiencies. Left column; yd — >■ rjnp 
(inclusive reaction), center column: 'yd — >■ rjp{n) (recoil protons 
detected), right column 7^ — >■ rjn{p) (recoil neutrons detected) 
for different incident photon energies as function of the r] 
cm polar angle. (Red) dots: efficiency from quasi-factorization 
(e/ac), (black) open squares: efficiency from phase-space simu- 
lation( Eph). 



Possible deviations might occur in the three-body final 
state of meson, participant, and spectator nucleon for ex- 
ample due to final state interactions. 

In the second approach, the above discussed distribu- 
tions of the detection efiiciencies for mesons and recoil 
nucleons as function of particle laboratory angle and ki- 
netic energy were applied event-by-event to the data. This 
step is completely model independent since the correction 
is only based on measured quantities (recoil nucleon ener- 
gies from time-of-flight in TAPS and from reaction kine- 
matics in CB). However, it does not correct for the missing 
mass cuts (roughly a factor of two for the cut at Am < 
MeV). The loss factor due to this cut was again simulated 
with an event generator based on reaction phase space and 
corrected. The results for this almost model independent 



efficiencies tfac from the quasi-factorization are shown as 
red dots in Fig. [121 Obviously, the results obtained with 
both methods are in quite good agreement. 

As expected, due to the trigger efficiency, for all reac- 
tion channels the detection efficiency is small at backward 
angles, in particular at higher incident photon energies. 
For the exclusive reactions, in particular the proton chan- 
nel, it is also small for meson forward angles, correspond- 
ing to recoil nucleons emitted at large angles with small 
kinetic laboratory energies. 




1000 1500 2000 ^ [UeM] ^°°° ''^°° ^°°° 



Fig. 13. Angle integrated total detection efficiencies eph (open 
squares) and e/ac as function of incident photon energy. 



The angle integrated detection efficiencies as function 
of incident photon energy are shown in Fig. 1131 The curve 
for the proton channel shows some structure at photon en- 
ergies of 1 GeV. This can be explained with Figs.[7l[T0l At 
the lowest incident photon energies, all recoil protons are 
detected in TAPS. Around 900 MeV, part of the protons 
reaches the Barrel. Since the detection threshold in the 
Barrel is higher than in TAPS, this leads to a decrease 
of the efficiency. But since proton detection above the 
threshold was more efficient in the Barrel («90%) than in 
TAPS (Ri60%) the overall efficiency rises again for higher 
incident photon energies. Since the neutron detection ef- 
ficiency varies smoothly for both detectors no structure is 
present for quasi-free production off the neutron. 



4.3.3 Summary of systematic uncertainties 

Typical statistical uncertainties for the total cross sections 
range from 0.5 % - 5 % for the inclusive data, from 1 % - 
10 % for the data with coincident protons and from 2 % - 
20% for the data with coincident neutrons (the first num- 
ber corresponds to the range of the Sii(1535) peak max- 
imum, the second number to maximum incident photon 
energies). As discussed below, systematic uncertainties are 
of comparable size for high incident photon energies but 
dominate in the Sn range. 

Systematic uncertainties are in three different cate- 
gories: overall uncertainties which cancel exactly in the 
comparison of different reaction channels, uncertainties 
which are similar for different reaction channels and can- 
cel to a large extent in ratios, and reaction channel related 
uncertainties which do not cancel. 

In the first category all cross sections are subject to 
an overall systematic uncertainty of the photon flux of 
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«10 %. The overall uncertainty of the target thickness of 
a few per cent is comparably small and the uncertainty of 
the decay branching ratio (below 1 %) is negligible. 

The second class of uncertainties is related to the de- 
tection of the ?7-mesons. The main steps of the analysis, 
which have to be reproduced by the Monte Carlo simula- 
tions, are the detection, identification, and calibration of 
photon showers, the invariant mass analysis for the iden- 
tification of the mesons, and the missing mass analysis re- 
moving background from reactions with additional mesons 
in the final state. 

At the most basic level, a stringent limit for system- 
atic uncertainties arising from the detection of photons 
and the subsequent identification of mesons via the in- 
variant mass analysis can be derived from a comparison 
of the results for 77-photoproduction off the free proton. 
This has been previously analyzed for the same setup for 
the ?7 — J> 27 and the ry — >■ 67 decay channels P51[^B5] . 
Systematic effects would enter cubed into the six-photon 
channel, but agreement was found on average at the 2 % 
level. Furthermore, also the present simulations with the 
two different methods - one relying on the phase-space 
event generator, the other using the quasi-factorization of 
the detection efficiency, are in good agreement. From this 
we estimate a typical 5 % uncertainty for a successful 77- 
reconstruction. The uncertainty related to the separation 
of signal and background in the invariant mass spectra 
(fitted with line shape and background polynomial) is not 
explicitly treated as an additional systematic effect, but 
the fit uncertainty is included into the statistical errors. 

Finally, the missing mass cut deserves special attention 
(see Fig. in]). At incident photon energies below 0.8 GeV, 
the spectra are practically background free. At higher en- 
ergies background from rjir final states and the tail of the 
missing mass distribution for single ry-production arising 
from the momentum distribution of the bound nucleons 
mix. The simulations indicate that the background reac- 
tions do not contribute in the region of negative missing 
masses (cf. Fig.IB]). Since only those events were accepted, 
background contamination is estimated at most at the per 
cent level. However, due to this cut, the simulation must 
closely reproduce the shape of the missing mass peak, in- 
cluding effects of Fermi motion. This is the case for the 
background free peaks at low incident photon energies and 
for all energies for the peak shape at the non-contaminated 
side. In the background contaminated regions the data can 
be reproduced by a summation of the simulated structures 
for peak and background (with properly adjusted relative 
contributions), but no direct check is possible. From the 
agreement between data and simulations and the fraction 
of the missing mass signal extending into the tail region 
we estimate for the missing mass analysis a systematic un- 
certainty of 3 % for the total cross section in the threshold 
region rising to 15 % at 2.5 GeV. Altogether, independent 
on the reaction channel we estimate the systematic uncer- 
tainty for the 77 detection at 8 % at production threshold 
rising to 20 % at 2.5 GeV. 

The most critical uncertainties are related to the detec- 
tion of the recoil nucleons. From the agreement between 



Table 2. Summary of systematic uncertainties for the quasi- 
free reactions. ^' photon flux, target thickness, decay branching 
ratios, ^' trigger efficiency, 77 analysis cuts, 77 detection effi- 
ciency. When two numbers are given the first corresponds to 
the threshold energy range, the second to E^— 2.5 GeV, and 
linear interpolation is used in between. 



source 


7d -^ {n)pT] 


"fd — S- (p)7l77 


overall normalization^-' 


10% 


10 % 


77 detection^-' 


8 % - 20 % 


8 % - 20 % 


recoil nucleon detection 


8% 


15 % 



simulated and measured efficiencies we estimate on aver- 
age ~8 % uncertainty for recoil protons and 15 % for recoil 
neutrons. Uncertainties can be larger (in particular for the 
proton) for kinematical parameters where they are partly 
detected in TAPS and partly in the CB, i.e. for incident 
photon energies from 0.8 GeV - 0.9 GeV and cos(6'*)> 
and at higher incident photon energies for cos(6'*)?s -0.1 
(cf. Fig. [T]) . The above discussed systematic uncertainties 
are summarized in Tab. [5J 
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Fig. 14. Comparison of total cross sections. (Blue) upward tri- 
angles: quasi-free proton cross section ap, (red) dots: quasi-free 
neutron cross section cr„, (black) open squares: inclusive quasi- 
free cross section a„p, (black) stars: (t„ + ap. Downward (ma- 
genta) triangles: inclusive quasi-free cross section from Weiss 
et al. '351. Insert: ratio of neutron cross sections. 



As already discussed in P^fTS] the nucleon detection 
uncertainty may be checked in an independent way. The 
cross section of the coherent process ■^d — )• ^77 is negligible 
PO] compared to the quasi-free reaction. Therefore, the 
quasi-free reactions cr„p (inclusive, no condition for recoil 
nucleons), ap (coincident recoil protons), and <7„ (coinci- 
dent recoil neutrons) must obey (T„p = <Jp + cr„. This is 
demonstrated in Fig. [TH where the sum of the quasi-free 
proton and neutron cross sections is compared to the in- 
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elusive cross section. The agreement is excellent and allows 
an independent extraction of the neutron cross section as 
a'^ = anp — <Jp- The ratio (Jn/cr'n of the two results is shown 
in the insert of Fig. [TH The agreement is within statistical 
uncertainties for most data points and typical deviations 
do not exceed the 15 % level. As a further test, the distri- 
bution of the deviations Sai normalized by the statistical 
uncertainties Aa^ 



Act, 



da'^/dfi — dunl dfi 

y/A'^{da'„/d[2) + A^{dan/df2) 



(9) 



for all data points (420 entries) of the angular distribu- 
tions from production threshold to 2.5 GeV is compared 
to a Gaussian distribution in Fig. [121 The fitted Gaussian 
distribution corresponds to a width of u = (1.25 ± 0.10) 
and a mean oi ^ — (0.034 ± 0.110), fairly close to a stan- 
dard Gaussian distribution. In particular, the mean is not 
significantly different from zero so that no indication for 
a systematic deviation is indicated. 
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Fig. 15. Distribution of deviations between dan/dO and 
da'n/dn = danp/dn — dap/dQ. Solid (red) curve: fitted Gaus- 
sian distribution (width a=1.25, mean /x=0.034), dashed (blue) 
curve: standard Gauss: ((t=1, /i=0) 



This is a very stringent test for systematic uncertain- 
ties related to the recoil nucleon detection since it is the 
neutron detection efficiency which enters in the extraction 
of (T„, while only the inherently different proton detection 
efficiency enters into tr^. The result means that the corre- 
sponding systematic uncertainties quoted in Tab[2]derived 
from the analysis of the recoil nucleon detection efficien- 
cies are probably pessimistic. 

Due to this good agreement total neutron cross sec- 
tions as function of incident photon energy are given in 
this paper as weighted averages ((t„) of Un and a'^ which 
improves the statistical quality. For the shape of the an- 
gular distributions only the direct measurement with the 
neutron coincidence is used since due to the small detec- 
tion efficiency at extreme angles (in particular for the pro- 
ton) systematic effects for a'^ are larger. 



5 Results and discussion 

Throughout this paper all quasi-free differential cross sec- 
tions are given in the cm (center-of-momentum) system of 
the incident photon and a target nucleon at rest. This sim- 
plifies the comparison to angular distributions measured 
off the free nucleon since apart from the immediate thresh- 
old region such quasi-free cross sections are only moder- 
ately smeared out by Fermi motion in this system, while 
in the ^d system they have completely different shapes 
due to the Lorentz boosts (see ^36^ for details). 



5.1 Comparison to previous results 

The inclusive cross section cr„p has been previously mea- 
sured twice [36 , 39 at the Mainz MAMI accelerator with 
different configurations of the TAPS detector (the setup 
used for ^39^ covered a larger solid angle than for [3B] ) up 
to incident photon energies of 0.8 GeV. The total cross 
section from [39] is included in Fig. [141 Typical angular 
distributions are compared in Fig. [T51 
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Fig. 16. Comparison of the inclusive cross section (j„p to pre- 
vious results (Krusche95 [36], WeissOS [39]). Errors are only 
statistical. 



The shapes of the angular distributions are similar for 
all three experiments. There is a systematic deviation be- 
tween the two previous measurements with the present 
data for the absolute scale of the two lowest energy bins. 
These are, however, probably not due to the normalization 
of the cross section data but to the systematic uncertainty 
in the measurement of the incident photon energy. Both 
bins are located in the extremely sharp rise of the cross 
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section close to threshold (of. Fig. [H)) . In this range, al- 
ready small effects in the determination of the incident 
photon beam energy are strongly amplified in the magni- 
tude of the cross section. The resolution for the incident 
photon energy in the threshold region was better for the 
previous MAMI experiments, which aimed at very precise 
threshold measurements in view of FSI effects. The MAMI 
data had typically 2 MeV bin width for the incident pho- 
ton energy with an absolute calibration uncertainty of less 
than 1 MeV while the present data were measured with 
10 MeV bin width and a calibration uncertainty of several 
MeV. Therefore, the systematic quality of the previous 
data in the immediate threshold region is almost certainly 
superior. However, this region is not of much interest for 
the present experiment. Agreement is much better and 
within systematic uncertainties at higher incident photon 
energies. In this range, the data from |3S] have somewhat 
larger systematic uncertainties than the other two data 
sets due to the restricted solid angle coverage. 



of free and folded cross section (solid line in the insert 
of Fig. ITTI) were then applied as correction factors to the 
measured quasi-free cross section. Measured and corrected 
quasi-free cross sections are compared in the insert of Fig. 
[T71 The corrected quasi-free cross section is then com- 
pared in Fig. [T7] to the data base of free proton results. 
The agreement between the free and quasi-free data is ex- 
cellent for most incident photon energies. This indicates 
that no nuclear effects (FSI, re-scattering) except nuclear 
Fermi motion influence the quasi-free data. The latter ef- 
fects are significant in the steep slopes of the cross section 
for photon energies below 1.1 GeV, they are negligible in 
the flat region at higher incident photon energies. 

Angular distributions for the free and quasi-free 7p ^• 
pi] reactions are summarized in Fig. 1181 The quasi-free 
data have not been corrected for Fermi motion. They have 
been fitted with Legendre polynomials 
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Fig. 17. Comparison of the total cross section of quasi-free 
production off the bound proton to free proton data. Filled 
(red) circles: quasi-free (jp corrected for effects of Fermi mo- 
tion (see text). Open (blue) diamonds: free proton data from 
[roi[2T1E2 , 23, 28, 2 51151] . Total cross sections for the free proton 
data from 22,29 have been estimated from the published dif- 
ferential cross sections. Solid line: eta-Maid model [53]. Insert: 
comparison of quasi-free proton data (black triangles) to quasi- 
free proton data after correction of Fermi motion effects (red 
filled circles) (see text). Solid curve: ratio of free and folded 
cross section (scale at right hand side). 



The quasi- free proton data are compared in Figs. 11711^ 
to free proton data. The total cross section for the quasi- 
free reaction has been corrected for the effects of Fermi 
motion in the following way. The well known energy de- 
pendence of the cross section for ?7-production off the free 
proton was folded with the momentum distribution of the 
bound proton using the deuteron wave function in mo- 
mentum space from [68| as described in [36) . The ratios 



|^-|EA«(».(e;)), 



(10) 



where the Ai are expansion coefflcients (higher order coef- 
ficients were not significant). The phase-space factor qt/k* 
(<Zn, k*'- meson and photon cm momenta) is evaluated for 
the photon - nucleon-at-rest cm system. 

At incident photon energies above 1 . 1 GeV the angular 
distributions are in excellent agreement with the quasi-free 
data. Close to threshold the large influence of Fermi mo- 
tion is visible, but the comparison of the quasi-free data to 
the MAID model result folded with Fermi motion (dashed 
lines) demonstrates that this effect is well under control. 
A large deviation between free and quasi-free data occurs 
also for incident photon energies around 975 MeV. In this 
region the pronounced 'dip' in the total free cross section 
is fllled in by Fermi motion from the tail of the Sn res- 
onance, however again the folded cross section is in good 
agreement with experiment. 

In summary, we conclude that after correction for the 
effects of nuclear Fermi motion the absolute scale and the 
shape of the angular distributions of the quasi-free proton 
data agree very well with the most recent and most precise 
measurements of 77-photoproduction from the free proton. 
This is the systematic basis for the discussion of the quasi- 
free neutron data. 



5.2 The quasi-free reaction •yn — >■ nrj off the neutron 

The total cross section of the quasi-free reaction off the 
neutron is compared to the quasi-free proton data in Fig. 
[T9I which also shows the neutron/proton cross section ra- 
tio. The behavior at low incident photon energies {E^ <800 
MeV) is consistent with previous results. The cross section 
ratio in the Sn-region is close to 2/3 and rises to the kine- 
matic threshold close to unity because in the immediate 
vicinity of the threshold the participant - spectator ap- 
proach becomes meaningless (dictated by energy and mo- 
mentum conservation at threshold 'participant' and 'spec- 
tator' nucleon have always identical momenta). 
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Fig. 18. Differential cross sections for 7p — > prj. (Red) filled circles: quasi-free data from present experiment (not corrected 
for Fermi motion), (Red) solid lines: fits with Legendre polynomials. Other symbols free proton data: (blue) squares: [28) . 
(black) stars: [29) . (green) triangles: [31] (data have been partly re-binned to cover the same energy ranges). Dashed (blue) lines: 
Eta-Maid model [33] folded with Fermi motion 
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Fig. 19. Comparison of quasi-free proton and neutron excita- 
tion function. Upper part: Curves model results, dashed: Eta- 
MAID for proton ^33j, solid: Eta-MAID for neutron ]33j, dash- 
dotted: Shklyar et al. [49]. Insert: comparison of the total neu- 
tron cross sections extracted from the coincident measurement 
of neutrons (un) and the difference of inclusive and proton data 
(a'n)- Bottom part: Cross section ratio a„/ap compared to pre- 
vious data from quasi-free production off the deuteron [32 and 
off *He [38] and model results (solid: MAID [33], dashed: Shkl- 
yar et al. 33]) folded with Fermi motion. 



Around incident photon energies of 1 GeV, correspond- 
ing to W « 1.7 GeV a bump-like structure is visible in the 
neutron cross section, which is not seen for the proton. 



In fact, it is even more pronounced in the ratio of neu- 
tron and proton data, which shows a sharp rise around 
Ey—1 GeV. Before we discuss this structure in detail, we 
compare the angular distributions of the two reactions, 
which are summarized in Fig. 1201 They have been fitted 
with the Legendre series from Eq. [TU] and the coefficients 
are shown in Fig. [21] 

The coefficients of the quasi-free proton data are in 
good agreement with the free proton data as expected 
from Fig. [18] The largest deviations occur for a few values 
of the A3 coefficient close to the threshold, this might be 
due to uncorrected Fermi motion effects. 

The neutron results obtained from the analysis with 
coincident recoil neutrons ((T„) and from the difference 
of inclusive data and data with coincident protons cr,^ = 
Cnp — cTp (filled and open circles in Fig. 1201) are also in 
good overall agreement, some discrepancies occur for the 
extreme angles, in particular in forward direction. Here, 
one should note that the detection efficiency for the re- 
action with coincident recoil protons almost vanishes for 
the extreme forward angles of the 7y- meson (see Fig. [T2]) . 
so that a'^ is less well defined in this regime. 

For the comparison of proton and neutron cross sec- 
tions we discuss three different energy ranges. At high in- 
cident photon energies above 1.5 GeV, the absolute mag- 
nitude of the cross sections as well as the strongly for- 
ward peaked shape of the angular distributions are almost 
identical. This is the energy region, were previous model 
analyses of free proton data (see e.g. [53]) have identi- 
fied dominant contributions from t-channel background 
terms. These contributions seem to be similar for protons 
and neutrons, which is not unexpected when considering 
isospin invariance. Agreement with the MAID model at 
higher incident photon energies is not good, the observed 
forward peaking of the cross section is not reproduced. 
The analysis of the previous ELSA proton data 23,251125] 
in the framework of the Bonn-Gatchina model had found 
a strong contribution of a Di5(2070) state to r/-production 
at large incident photon energies. This state is not visible 
as a bump in the total cross section, not even like the 
small indication of the P-resonances at W Kil.7 GeV). It 
was extracted from the analysis of the angular distribu- 
tions. The present neutron/proton ratio of the total cross 
section might show some structure in this energy region 
(cf. Fig. [ini bottom part), however at the very limit of 
statistical significance. 

The total cross section for quasi-free 77-production off 
the proton and the neutron are quite different at low inci- 
dent photon energies (below 900 MeV) in the range of 
the second resonance region. In this regime, the abso- 
lute magnitudes of the cross sections refiect the ratio of 
the helicity amplitudes of the Sii(1535) resonance (see 
Eas. lll2l) . The shape of the angular distributions is dom- 
inated by the interference between the Sii(1535) and the 
Di3(1520) resonances, which involves a term proportional 
to the A2 coefficient of the Legendre series. The electro- 
magnetic helicity-1/2 couplings A1/2 of the D13 have iden- 
tical signs for proton and neutron, while the Sn couplings 
have different signs. 
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Fig. 20. Quasi- free angular distributions, labels indicate incident photon energy. (Blue) squares: proton coincidence ap, (red) 
filled circles: neutron coincidence a„, (red) open circles: difference of inclusive and proton a'„. Dashed (blue) curves: fit of proton 
data, solid (red) curves: fit of neutron data, dotted (black) curves: Eta-MAID for neutron folded with Fermi motion. 
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Therefore, the interference term changes sign from proton 
to neutron, giving rise to negative A2 coefficients for the 
proton and positive ones for the neutron [3S]. This effect 
is also reflected in the model calculations [5511491150] . 

In the most interesting range around 1 GeV incident 
photon energy, where the peak-like structure appears in 
the total neutron cross section, proton and neutron an- 
gular distributions are not very different (apart from the 
absolute scale). In this region, the Sn - D13 interference is 
not visible any more {A2 for proton and neutron is small 
and negative) and the Ai coefficient shows a zero crossing 
with very steep rise. It was already discussed in |27 i that 
the simplest explanation for a rapidly varying Ai coeffi- 
cient is an interference between S and P waves. This ex- 
planation seems to be natural since it is well known that in 
this energy region the tail of the Sii(I535) resonance, the 
Sii(1650) resonance, and the Pii(1710) and/or Pi3(1730) 
resonances contribute. The model analyses still disagree 
in the relative importance of the two P-wave states. The 
77-MAID model [33] finds a dominant contribution from 
the Pii, while the Bonn-Gatchina analysis [S^ prefers 
the P13 state. Very preliminary results from a measure- 
ment of the 7p — > prj reaction (circularly polarized beam, 
longitudinally polarized target) at ELSA [5^ indicate a 
dominant helicity-1/2 component in this energy region, 
which is more in line with the MAID analysis favoring the 
Pii contribution. If one assumes only Sn and Pn states 
(£^0+ and Ml- multipoles), the Ai coefficient would sim- 
ply be proportional to Re(iJ*+ Mi- ) and thus the fast zero 
crossing would imply a rapid phase change between these 
multipoles indicating, that one is going through resonance 
around 1 GeV. 

It is, however, not yet understood, whether the struc- 
ture observed in the total neutron cross section at the 
same incident photon energy is somehow related to con- 
tributions from these resonances. Structures in the total 
cross section can obviously not arise from interferences 
between different partial waves. Therefore, different sce- 
narios involving contributions from specific resonances as 
well as interference patterns in the same partial wave have 
been discussed in the literature. Fix, Tiator, and Polyakov 
[52] have investigated whether the data could be consis- 
tent with the excitation of a narrow Pn-state. This work 
was motivated by the idea that the Pn-state of the pro- 
posed anti-decuplet of pentaquark states should be rel- 
atively narrow (width on the order of 10 MeV), have a 
strong electromagnetic coupling to the neutron, and a 
large rjN decay branching ratio !44,45]. They used two 
different versions of the 77-MAID model, the standard ver- 
sion [33] and the reggeized version [7D] as basis of their 
fits. The standard version, including a large contribution 
from the Di5(1675) resonance, reproduces fairly well the 
experimental ratio of neutron and proton cross sections, 
although it does not show the structure observed in the 
neutron data around photon energies of 900 MeV (see also 
discussion in the next section). The reggeized version with 
a much smaller contribution of the D15 reproduces the 
data only when an additional narrow resonance is intro- 
duced (taken as Pn). Since the data are smeared out by 
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Fig. 21. Coefficients of the Legendre series of the fits in 
Figs. I18I20I Left hand side: quasi-free proton. Open (red) cir- 
cles: present quasi-free data, (blue) downward triangles: Bonn 
ELSA data [21], (black) stars: CLAS data 29_, (green) up- 
ward triangles: Mainz MAMI data [31]. Right hand side: quasi- 
free neutron (from neutron coincidence (t„). Open (red) circles: 
present quasi- free data, solid lines: Eta- MAID [33] . Note the 
logarithmic scale for Aq,. 



Fermi motion, the width of this additional state is uncer- 
tain and could be as narrow as 10 MeV (roughly 40 MeV 
as upper limit) [52] . Similarly, an analysis performed in 
the framework of the Bonn-Gatchina model (BoGa) [53] 
can reproduce the neutron data reasonably well with three 
completely different scenarios, by either adding a 'conven- 
tionally' broad Pn resonance, a very narrow Pn state, or 
even by a careful adjustment of the interference pattern 
for the 5'- wave amplitudes. Shklyar, Lenske, and Mosel 
[in] find solutions in the framework of the Giessen cou- 
pled channel model with bump-like structures in the neu- 
tron excitation function around E~^ « 1 GeV just from 
coupled channel effects in the Sn - Pn sector, without in- 
troducing any additional resonance. A similar result from 
a coupled-channels K-matrix approach was presented by 
Shyam and Scholten [50] who report a bump-like struc- 
ture arising from superpositions and interferences of con- 
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tributions from the Sii(1535), Sii(1650), Pii(1710), and 
Pi3(1720) states. However, this structure appears broader 
than our experimental results discussed below. Finally, 
Doring and Nakayama |51j , using an 5- wave coupled chan- 
nel model, find a 'dip-bump' structure in the neutron cross 
section related to the opening strangeness thresholds of 
KA and KE photoproduction around 900 MeV and 1050 
MeV. Such unitary cusps are for example well-known for 
pion production reactions. The cusp structure in tt*' and 
TT"*" photoproduction off the proton at the 77-production 
threshold was discussed in detail in [7T|, and the cusp 
structure in jp -^ pn^ at the jp — > ?i7r+ threshold was 
analyzed in |72j. 

From the experimental side obviously two pieces of in- 
formation are missing to distinguish between these differ- 
ent scenarios. In the absence of any results for polariza- 
tion observables, it is impossible to isolate the responsible 
partial wave(s). Measurements of polarization observables 
like the helicity asymmetry E (longitudinally polarized 
target, circularly polarized beam), the target asymmetry 
T (transversely polarized target), and the asymmetry F 
(transversely polarized target, circularly polarized beam) 
have been initiated at the Bonn and Mainz accelerators, 
for 77 photoproduction off the free proton also at the CLAS 
experiment at Jlab. Furthermore, the influence of the mo- 
mentum distribution of the bound neutron obscures the 
intrinsic shape of the bump structure for the free neutron. 
To overcome this difficulty, we present in the following sec- 
tion a new analysis of the present data, based not on the 
incident photon energy but on the invariant mass of the 
rj — N final state, which is not affected by Fermi motion. 



5.3 The 77 - nucleon invariant mass distributions 

In the previous sections we discussed the cross sections 
as function of the incident photon energy measured with 
the tagging spectrometer. Due to the momentum distri- 
bution of the bound nucleons, each value of incident pho- 
ton energy corresponds to a broad distribution of invari- 
ant masses W of the rj - participant-nucleon pairs, giv- 
ing rise to the Fermi smearing of all narrow structures. 
However, in principle we can directly extract W from the 
four-momenta of the 77 and the participant nucleon. We 
have already shown in [Tl] that the bump around 1.7 GeV 
in W becomes then much more narrow. In that analysis 
only data with 77-mesons emitted at backward angles were 
included, since then the kinetic energy of the neutron de- 
tected in the TAPS detector can be determined from a 
time-of-flight measurement. However, this analysis can be 
extended to the full data set, using the kinematical overde- 
termination of the data. 

All kinematic variables (incident photon beam energy, 
target deuteron at rest) of the initial state are known. For 
the final state the four-momentum of the 7;-meson, the 
mass of participant and spectator nucleon and the recoil 
direction of the participant nucleon (polar angle O and 
azimuthal angle <?) are known. Missing is only the three- 
momentum of the spectator- and the kinetic energy of the 
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Fig. 22. Momentum distributions of spectator nucleons. (Red) 
dots: reconstructed from data, (black) lines: expected from 
deuteron wave function ^68!, (blue) squares: Monte Carlo sim- 
ulation including detector response. Left hand side: neutron 
spectator (i.e. recoil proton detected), right hand side: proton 
spectator (i.e. recoil neutron detected). 



participant nucleon. But these four variables can be ex- 
tracted using energy and momentum conservation which 
provide four equations. Monte Carlo simulations using the 
GEANT package have shown that with this reconstruction 
a typical experimental resolution of FWHM«25 MeV for 
W is achieved. 

As result of such an analysis Fig. [52] shows the dis- 
tribution of the momenta of the spectator nucleons, con- 
structed event-by-event from the reaction kinematics. In 
plane wave approximation with negligible FSI effects, these 
momenta must reflect the momentum distribution of the 
bound nucleons. As demonstrated in the figure, this is 
quite well fulfilled. Data generated with a participant - 
spectator Monte Carlo simulation using the deuteron wave 
function as input and including the response of the de- 
tector reproduce the measured distributions. For further 
analysis one can in principle cut away events with large 
spectator momenta, which are not close to quasi-free kine- 
matics. However, as it turned out (see Figs. [531 [M]), such 
a cut does not have much impact (in particular not on 
the angular distributions) apart from reducing counting 
statistics, so that it was not used for the differential cross 
sections. 

For the absolute normalization of the total cross sec- 
tion as function of W, the photon flux dN-y/dE-y measured 
as function of incident photon energy Ej , was folded with 
the nucleon momentum distribution to obtain the flux 
dN-y/dW in dependence of the final state invariant mass. 
The results for the total quasi-free cross sections of the 
proton and the neutron with and without cut on specta- 
tor momentum are summarized in Figs. [521 HH For both 
analyses good agreement between the quasi-free proton 
data and the world data set of free proton data is found. 
The neutron data show a pronounced, narrow peak around 
W fn 1.7 GeV. The position of this peak coincides with a 
dip in the proton excitation function. 
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In the following we discuss in more detail first the range 
of the Sii(1535) resonance and then the narrow structure 
in the neutron excitation function. 



5.3.1 The region of the Sii(1535) resonance peak 

The region of the Sii(1535) is interesting for two rea- 
sons. First since this is a well studied state, it can serve 
as a test case for the extraction of resonance parameters 
from quasi-free data with the above discussed kinemat- 
ical reconstruction. Furthermore, the results contribute 
to the discussion up to which energy range the observed 
cross section for 77-photoproduction is dominated by the 
Sii(1535). 
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where Wr, -Tr, are resonance position and width, fc^, g*^, 
q*jl are incident photon momentum, 77-momentum and tt- 
momentum at resonance position in the photon - nucleon- 
at-rest cm-frame, A^^ is the electromagnetic helicity-1/2 
coupling, bri—0.5, 67r=0.4, and 6^7r=0.1 have been used as 
branching ratios for the Nij, Nn, and Ntttt decays of the 
resonance, and rriN is the nucleon mass. 



Table 3. Result of Breit-Wigner fits. Nrj branching ratio of 
Sii is assumed as 6,,=0.5. Upper part of table: Comparison of 
fits of Sii(1535) resonance for free proton, quasi-free proton, 
and quasi-free neutron data to PDG estimates [34] and BoGa 
model fit [53] . f- Breit-Wigner mass, in brackets pole position; 
p. only magnitudes, no signs; §: only pole positions given, no 
Breit-Wigner mass. Bottom part: fit of neutron data with Sn. 
resonance and two further Breit-Wigner curves. All uncertain- 
ties of fit parameters are statistical only. 
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Fig. 23. Total cross sections as function of final state in- 
variant mass W without cut on spectator momentum. (Red) 
dots: quasi-free neutron, (blue) squares: quasi-free proton, 
(green) stars: free proton data. Curves: fitted (up to W = 
1600 MeV) Sii(1535) line shapes. (Black) solid: free proton, 
(blue) dashed: quasi-free proton, (red) dotted: quasi-free neu- 
tron. Dash-dotted curves: model results from 51 . Insert: ratio 
of quasi-free neutron - proton data. 



In the region of the Sii(1535) peak all data have been 
fitted with a parameterization of this resonance as Breit- 
Wigner curve with energy dependent width |12j : 
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The data have been fitted up to W Ril.6 GeV, where 
the proton and neutron cross sections start to deviate and 
obviously the line-shape can no longer be dominated by 
the Sii-resonance. The fit results are shown in Fig.E51and 
the parameters are summarized in the upper part of ta- 
ble|31 The results demonstrate the following. Breit-Wigner 
mass and width of the resonance extracted from the quasi- 
free proton and neutron data are in excellent agreement. 
The agreement with the free proton data is good, but not 
within statistical uncertainties. This had to be expected 
due to the finite W resolution of the quasi- free data, which 
tends to increase the width and to shift the resonance po- 
sition slightly upward. These parameters are also in good 
agreement with the values given by the particle data group 
[51] and the Bonn-Gatchina analysis [13] of the present 
data. Note that for the BoGa analysis not Breit-Wigner 
masses but pole positions are given, which agree with the 
PDG parameters. The almost perfect agreement of the 
shape of the Sn-peaks for the proton and neutron data 
in this range is also further evidence that this shape is 
alone dominated by the Sii(1535). Effects e.g. from the 
destructive interference of the two Sn-resonances, which 
are important at higher incident photon energies, should 
be different for proton and neutron since the ratios of 
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the electromagnetic couphngs of these two resonances are 
quite different for protons and neutrons. 

The electromagnetic helicity couplings A^,^ found from 
the fits agree for free and quasi-free proton data. The pro- 
ton couplings are slightly higher than the PDG value and 
the BoGa result. This could be a systematic effect, since all 
non-Sii(1535) contributions are neglected in the fit. The 
only large discrepancy arises for the neutron helicity cou- 
pling between the present and BoGa results on one hand 
side and the PDG value on the other side, which is much 
lower. Here, one should note that as already discussed in 
[12] there is a systematic discrepancy between the helic- 
ity couplings of the Sii(1535) extracted from pion pho- 
toproduction versus those from 77-photoproduction. The 
latter ones are significantly larger. The Sii(1535) dom- 
inates T] production but contributes only weakly to pion 
production which is dominated in this energy range by the 
Di3(1520). Therefore, 77-production is the better suited 
channel for the study of the Sii(1535) properties. In the 
meantime the PDG proton coupling became dominated 
by the larger values from 77 production, but the neutron 
coupling is still dominated by the small values from pion 
production. The resulting PDG neutron/proton ratio of 
the helicity couplings would correspond to a cross section 
ratio for rj production in the Sn maximum of 0.26, which 
is unrealistic. The BoGa analysis finds a ratio of 0.79 and 
the simple BW-fits a ratio of 0.68. 
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Fig. 24. Total cross sections as function of final state invari- 
ant mass W for spectator momenta ps < 100 MeV. Notation 
as in Fig. 1231 All curves for neutron data; dashed: fitted Sn 
line shape, dotted: broad Breit-Wigner resonance, dash-dotted: 
narrow Breit-Wigner, solid: sum of all. 



5.3.2 The region of the narrow peak in the jn -^ nrj 
reaction 

In order to estimate the width of the narrow structure 
observed in the neutron data, the excitation function has 
been fitted up to M^ w 1.8 GeV with a purely phenomeno- 
logical fit function. It is composed of the Breit-Wigner 
curve with energy dependent width for the Sii(1535) res- 
onance and two further simple Breit-Wigner curves with 
constant width {x =1). The curves are compared to the 
data in Fig. [M] and the fit parameters are listed in the 
bottom part of Tab. [31 The parameters obtained for the 
Sii are consistent with the results discussed above. The 
broad BW-curve located at M^ « 1.7 GeV just serves for 
the effective parameterization of the excitation function. 
It subsumes contributions from all normally broad reso- 
nances in this energy region (such as Pii(1710), Pi3(1720), 
D15I65O),...) as well as background components. The nar- 
row Breit-Wigner curve at M^ sal. 66 GeV has a FWHM of 
only (25±12) MeV, on the same order as the experimental 
resolution of 25 MeV (FWHM). This width is somewhat 
dependent on the chosen parameterization, but also trials 
with different background shapes, e.g. of polynomial type, 
which result in a poorer fit quality, indicate a width below 
the 50 MeV level. 

Since so far there is no information about the quan- 
tum numbers of this structure, in fact it is not even clear, 
whether it corresponds to a nucleon resonance, parame- 
ters like electromagnetic couplings cannot be given. How- 
ever, if we treat the structure as a narrow Sn resonance 



the normalization of the fit corresponds to A1/2 ■ \f^ ~ 
12 X 10-3GeV-i/2. 

The results for the angular dependence of the exci- 
tation functions are summarized in Figs. [HI [251 Due to 
statistical limitations in the extraction process of the cross 
sections depending on the final state W ^ the angular dis- 
tributions are only coarsely binned. Fig. [511 shows in the 
upper part excitation functions in dependence on W for 
four different bins of cm polar angles as well as the neu- 
tron/proton ratios. The bottom part shows the correspond- 
ing angular distributions for different bins of W . Finally, 
Fig. [2H1 summarizes coefficients of the Legendre series of 
Eq. [TOl fitted to the angular distributions. The results 
arc compared to free proton data and model calculations. 
The comparison to free proton data from the recent most 
precise measurement at MAMI [31] demonstrates impres- 
sively how well the elementary reaction on the free pro- 
ton can be approximated by quasi-free data with W re- 
constructed from the jtr] final state kinematics. Significant 
deviations occur only very close to the production thresh- 
old, where the effects from Fermi motion are most pro- 
nounced. Comparing proton and neutron data, in partic- 
ular the excitation functions for forward angles, it is even 
more apparent than in the total cross section that the nar- 
row structure observed in the neutron excitation function 
is accompanied by a pronounced dip in the proton data 
at the same position and of comparable width. It seems 
to be highly unlikely that these two structures are unre- 
lated. This might indicate that some interference with a 
sign change between proton and neutron is involved. 
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Fig. 25. Top, first row: excitation functions for different bins of rj cm polar angle. (Blue) open squares: quasi-free proton data, 
(black) stars: free proton data from [31], (red) dots: quasi-free neutron data scaled up by 3/2. (Blue) solid lines: 77- MAID [33j for 
the proton target, (red) dashed lines: r;-MAID for the neutron target. Second row: ratio of neutron and proton cross section for 
data and ?7-MAID. Vertical dotted lines: position of narrow peak in neutron data, horizontal dotted lines: a„/ap=2/3. Bottom 
part: angular distributions. Same notation as top part (but neutron data not scaled); additional curves (blue) dotted: proton 
model from Shklyar et. al. [35], (red) dash-dotted: neutron model from same reference. 
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Fig. 26. Fitted coefficients of tiie Legendre series as function 
of final state invariant mass W. Left hand side: (black) stars 
free proton, (blue) squares quasi-free proton. Right hand side: 
(red) dots quasi-free neutron. Curves: dashed standard MAID 
model [33], solid: Shklyar et al. [49] . 



However, since the neutron peak and the proton dip are 
also visible in the total cross section, at least part of this 
interference must be in the same partial wave. A dip-like 
structure has also been observed for comparable values of 
W in the 7r~p — > r]n reaction, although at much lower 
statistical significance. As a possible explanation, similar 
to one version of the BoGa-model [S3], the interference 
between the two Sn resonances was discussed in a K- 
matrix approach [75] . 

The comparison of the neutron data to model predic- 
tions leads to the following conclusions. Models which try 
to explain the structure observed in the Fermi smeared 
excitation function of the neutron data by one conven- 
tionally broad nucleon resonance like one of the scenarios 
in [S3] are ruled out by the narrow width (see Fig. [Ml) on 
the order of 25 MeV. 



The coupled channel approach of Doring and Nakayama 
[ST] shown in Fig. [23] in fact produces a structure of simi- 
lar width close to the observed position, although the ex- 
act shape is somewhat different. Since this model includes 
only s-wave contributions it cannot predict realistic shapes 
of angular distributions. Nevertheless, the fitted Legendre 
coefficients (see Fig. [26t do at least not contradict the as- 
sumption that only s-waves are important. The peak-like 
structure is clearly seen only in Aq, which is proportional 
to the total cross section. From the higher coefficients only 
A3 might show a little structure close by, albeit not sta- 
tistically significant. 

A direct comparison of the predictions of the MAID 
model [33] and the model from Shklyar et al. [JU] to the an- 
gular distributions (see bottom part of Fig. [25]) shows that 
both models disagree with the neutron data throughout 
most of the energy range, in particular around the peak 
structure at W ^ 1.67 GeV. Agreement with the proton 
data is of course much better since free proton data have 
been used to fix the parameters of both models. 

A more detailed comparison reveals some interesting 
features. The MAID model reproduces reasonably well the 
prominent structure in the ratio of neutron and proton 
cross section in the forward angular range (upper part of 
Fig. [25]). However, even there it does not show any indica- 
tion of the peak- like structure in the excitation function of 
the jn -^ nrj reaction. The peak in the ratio stems alone 
from the dip in the proton cross section. Since it is un- 
likely that these two structures are unrelated, this casts 
also some doubts whether the 'dip' structure in the proton 
cross section was correctly interpreted in this model. The 
comparison certainly rules out that the neutron structure 
can be entirely explained by a strong contribution of the 
Di5(1675) resonance. A comparison of the fitted Legendre 
coefficients shows that the MAID model does not repro- 
duce the peak in the Aq coefficients and also fails com- 
pletely for Al, indicating that the strong S—P interference 
is not reflected in the model. On the other hand, the good 
agreement for the A2 coefficient means that the S — D 
interference term (in particular the Sii(1535) - Di3(1520) 
interference) is well understood. The Shklyar model [?5] 
shows at least some indication for a peak in Aq at the 
right position and is in much better agreement with Ai, 
although here it predicts a dip structure around W «1.67 
GeV, which is not in the data. Agreement with A2 is not 
as good as for MAID. It is evident that the comparison of 
data and models does not allow a final conclusion about 
the nature of the structure in the neutron excitation func- 
tion. 



6 Summary and Conclusions 

Precise cross section data have been measured for quasi- 
free photoproduction of 77-mesons off protons and neutrons 
bound in the deuteron. Due to the combined analysis 
of events in coincidence with recoil protons, recoil neu- 
trons, and of the inclusive reaction systematic uncertain- 
ties related to the detection of the recoil nucleons could 
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be reliably controlled. The results confirm earlier mea- 
surements in the region of the Sii(1535) resonance and 
reveal an unexpected structure in the total cross section 
of the 7n -^ nrj reaction around incident photon energies 
of 1 GeV. 

The results of an analysis based on the invariant mass 
of the ?7-proton pairs from quasi-free photoproduction off 
the proton are in excellent agreement with free proton 
data. This demonstrates that no nuclear effects beyond 
Fermi motion are involved and that these effects can be 
reliably removed by this method. An identical analysis of 
the quasi-free neutron data confirms a very narrow peak 
in the neutron excitation function around W f«1.67 GeV 
with a width of only ^25 MeV (FWHM), which seems to 
correspond to a dip-like structure in the proton excitation 
function at the same energy. The nature of this structures 
is not yet understood. Clearly ruled out are only single 
isolated resonances with conventional width like e.g. the 
Di5(1675) in the MAID model. Scenarios with a narrow 
resonance or different types of coupled channel effects can- 
not yet be discriminated with the available data. 
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